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The least action principle is exploited as a simulation tool to find the optimal dynamic path 
for spatially extended systems driven by a small noise. Applications are presented for thermally 
activated switching of a spatially-extended bistable system as well as the switching dynamics of 
magnetic thin films. The issue of nucleation versus propagation is discussed and the scaling for the 
number of nucleation events as a function of the terminal time and other material parameters is 
computed. 
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Thermal and other noises play a very important role 
in the behavior of many systems. No matter how small 
the noise is, it always has a non-trivial effect on long 
enough time intervals. Thermal noise has an especially 
important role for micro- or nano-devices since the en- 
ergy barrier between metastable states may come close to 
UbT in such systems and it may in fact limit the smallest 
size of such devices. One example is found in magnetic 
recording industry where the superparamagnetic limit, 
below which the data retention time becomes too short 
for commercial purposes, is now considered a serious lim- 
iting factor for the maximally achievable storage density 
jjj. In order to better control such processes, it is neces- 
sary to first obtain a detailed quantitative understanding 
of the effect of the noises in such systems. For bistable 
systems that are often used in storage and memory appli- 
cations, issues such as the mean switching time and the 
dynamics of the switching path have to be addressed. A 
major challenge is that most systems of interest are spa- 
tially extended and switching cannot be described using 
only a few degree of freedom. For instance, magnetiza- 
tion reversal in micron-sized magnetic films proceeds by 
nucleation and domain wall motion instead of coherent 
rotation^) . Predictions of reversal rates based stochastic 
coherent rotation models jj| (stochastic ordinary differ- 
ential equations) are off by orders of magnitude In 
these cases, one has to deal with the added complexity 
of stochastic partial differential equations. 

Traditionally the methods of choice for a quantitative 
understanding of the effect of noise have been the Monte 
Carlo method or direct simulation of the Langevin equa- 
tion. For small level of noise, however, these methods 
become prohibitively expensive. Wile the maximal time 
step for Monte Carlo or Langevin simulations is still re- 
stricted by the deterministic part of the dynamics, the 
effect of the noise is significant only on time scales which 
are exponentially large in the inverse of noise amplitude. 
This time scale is the natural one for the noisy dynam- 
ics, and it is unaccessible by Monte Carlo or Langevin 
simulations. For instance, for the magnetic problem (of 
interest in magnetic recording) , the times scale accessible 
to Monte Carlo or Langevin is on the order of nanosec- 



onds, and the time scale of interest might be years. 

In this paper we introduce a new numerical procedure 
to study spatially extended systems driven by a small 
noise. We are interested in the effect of rare events over 
long periods of time. Our method is based on the theory 
of large deviations which provides a least action princi- 
ple for the most probable dynamic path Q. Even tough 
there has already been works using the least action prin- 
ciple for ordinary differential equations (see e.g. |6|), ex- 
tending this in the context of spatially extended systems 
present non-trivial numerical difficulties. In the exam- 
ples we consider, namely thermally activated switching 
in a system with a bistable potential modeled by the 
Ginzburg-Landau equations in one and two dimensions 
and sub-micron-sized magnetic thin films, the numeri- 
cal computations are non-trivial because switching pro- 
ceeds by nucleation and domain wall motion which are 
both very localized in space and time, and require using 
very fine grids. Also, the diffusion terms present in these 
equations lead to very bad condition numbers, which re- 
sults in a loss of accuracy as well as an increase of the 
numbers of iterations necessary for convergence. We in- 
troduce an efficient numerical method to overcome these 
problems which incorporates amongst other techniques 
the use of suitable preconditioners, efficient minimiza- 
tions techniques such as the Quasi-Newton method, etc. 

Let us begin by considering the simplest continuum 
model describing the spatio-temporal states of a bistable 
system, namely the Ginzburg-Landau equation 

u t = Su xx -S- 1 V'(u), (1) 

on the interval [0, 1] with the boundary conditions u(0) = 
u(l) = 0. We take V to be a standard double- well poten- 
tial, V(u) = j(l — u 2 ) 2 . (|lj) is expressed in appropriate 
dimensionless variables in which 6 is a small parameter 
which indicates that the reaction term, S~ 1 V'(u), is fast 
while the diffusion is slow. (|l|) can be considered as the 
gradient flow associated with the energy 

E{u] = ±j\sul + 26- 1 V(u))dx. (2) 
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FIG. 1: Snapshots of profiles of the minimizer u during a 
switching from u+ (top curve) to U- (bottom curve) at six 
equally spaced times on [0, T] for different T at fixed S = 0.03: 
(a) T = 7; (b) T = 2; (c) T = 1; (d) T = 0.8; (e) T = 0.6; (f) 
T — 0.4. As T decreases, the number of nucleations increases: 
one in (a), two in (b), (c), and (d), three in (e) and (f). 

The dynamics in (Q) has two stable equilibrium states, 
u + and U-, which minimize the energy (^). When 5 is 
small, u±{x) = ±1 except at a thin boundary layer of 
width 5 at x = 0, 1 (see Fig. 0). 

Now let us add to (||) a small noise term modeling 
thermal effects: 

u t = Su xx - 5~ l V'(u) + yfsr], (3) 

where e is proportional to the temperature of the system 
and T] is space-time Gaussian white noise with covariance 

(r)(x,t)r](y,s)) = 5(x -y)S(t - s). (4) 

The presence of the noise in (||) destroys the long-time 
stability of the equilibria u±. For instance, if the initial 
state is u+, there is a finite probability that the system 
switches to u_ in any time interval [0, T]. To quantify 
the probability of such a switching or, in fact, any other 
event, one can introduce an action functional, Sr[u], such 
that for small e, the probability that the solution u of (||) 
be close to a given path ip on [0, T] can be estimated by 
§0 

Prob{uw (/ 3}-exp(-e~ 1 S , T M). (5) 

The action functional corresponding to (^) is given ex- 
plicitly by 

S T [u}= / (u t -5u xx + 5- 1 V'(u))dxdt. (6) 
Jo Jo 

The large deviation principle in (||) allows us to esti- 
mate the probability of various events associated with (0) 
by constrained minimization of the action functional m). 
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FIG. 2: Values of the action as a function of switching time 
T for the six minimizers shown in Fig. 1: the bottom curve 
corresponds to (a) in Fig. 1, while the top curve corresponds 
to (f ) . Each switching path it is a local minimizer of the action 
for all T: however, the glnobal minimizer changes with T. 

For instance, the probability Pt that the system switches 
from m + to u- before time T satisfies 

lim £ In Pt = — minjS'o tM} , (7) 

where the minimization in (Q) is constrained by the 
boundary conditions 

u|a!=o = «U=i = 0, u\ t=0 = u + , u\ t=T = u-. (8) 

Furthermore assume that the system switches before time 
T . Then the minimizer of (Q) gives the optimal switching 
path between u + and u_ during the time interval [0, T], 
in the sense that the probability that the system switches 
by another path is exponentially smaller in e. 

To find the optimal path, we numerically minimize the 
action functional, subject to the boundary condition in 
(^J). We put down a numerical grid on the space-time 
domain (0,1) x [0, T}. The action functional (^J) is dis- 
cretized using finite difference formulas. Numerical opti- 
mization is implemented using a Quasi-Newton method, 
the BFGS method ||. In order to speed up convergence, 
we used the operator obtained from minimizing the linear 
part of (|^), J T (u t —5u xx ) 2 dxdt, as the preconditioner. 
Similar procedures were used for the other examples be- 
low and numerical details are provided in [ p~0| | . 

Figure |l] shows the sequence of profiles of u at differ- 
ent times in [0, T] for various values of T at a fixed <5 
||. The switching proceeds by nucleation followed by 
propagation of domain walls. For large T, the switch- 
ing proceeds by propagating a domain wall from one 
boundary to the other. As T becomes smaller, how- 
ever, the number of nucleation events increase. Such 
switching scenario with multi-nucleations are relevant at 
finite temperature fiofl . In Figure ^, where we display 
the values the action for the various local minimizers, 
plotted against T. Figure || shows the space-time plot of 
(Ut — Su xx + S^ 1 V'(u)) during a switching event. This 
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FIG. 3: Space-time value of (u t - SAu + 5~ V'(ii)) 2 for the 
minimizer (f) in Fig. Q. This can be interpreted as the min- 
imal noise necessary for switching: the noise induces the nu- 
cleations then propagates the domain walls. The peaks cor- 
respond to nucleations. 



corresponds to the minimal (squared) noise necessary to 
make a switching. 

These results can be understood as follows. Consider 
the critical points of the energy (0), i.e. the solutions of 
= 8u xx — ($ _1 V'(u) with m| k= o = w| a= i = 0. Besides u+ 
and it_, corresponding to the minimizers of the energy, 
there are also saddle point configurations, denoted by 
us, with an increasing number of domain walls. For each 
saddle point, there is a path joining u + to u_ through 
that saddle point. Our result shows that, for large T, 
the switching path crosses the saddle point configuration 
with minimum energy, i.e the configuration with a single 
domain wall (see (a) in Figure [l]). As T is decreased, 
one might naively think that the switching will occur via 
the same path followed at a faster rate. This intuition is 
wrong. For smaller T, the optimal switching path crosses 
a saddle point configuration with increasing energy, i.e. 
more nucleations and therefore more domain walls, givin, 
rise to the cascade of nucleation events seen in Figure 
The reason is simply that both nucleation and domain 
wall motion are noise induced. As T decreases, at fixed 
number of nucleations the speed of propagation of the 
domain wall must increase in order to achieve complete 
switching during the allowed time T. This is energy con- 
suming and, at certain critical values of T it becomes 
more favorable to make an additional nucleation. 

The same type of cascade in the number of nucle- 
ations is also observed if T is kept fixed but 8 is de- 
creased, because the cost of domain- wall motion increases 
as 8 decreases, whereas the cost of a nucleation stays 
the same. In fact, a simple scaling argument gives 
n* = C(i5T) _1 / 2 , for small 8 and T, where n± is the 
number of nucleations in the optimal switching path, 
and C is a numerical constant. Furthermore, one has 
linwjelnPT = -C'iST)- 1 / 2 , which gives the envelop 
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FIG. 4: Snapshots of profiles of the minimizer u during a 
switching from u+ (top figure in one column) to u~ (bottom 
figure in one column) at five equally spaced times on [0, T] 
for T — 1 and different 8 and different boundary conditions: 

(a) -(c) correspond to case (i), (d) to case (ii). In (a) 8 = 0.04; 

(b) S = 0.03; (c) 8 = 0.01; (d) 8 = 0.01. The gray-scale is 
from white for u = —1 to black for u = +1. 



of the curves in Fig. |2[ The detailed analysis of the cas- 
cade process of nucleations will be presented elsewhere. 
As an example in two dimensions, we minimize 



St[u] 



{ut-8ku + 8- x V'{u)) 2 dPxdt, (9) 



where Vt is the unit square, il = [0, 1] x [0, 1]. We will 
present results with two different boundary conditions: 
(i) u = 1 at x = and x = 1, u — — 1 at y = and 
y = 1, and (ii) u = on the edge of the square. In both 
cases there are two global minimizers of the energy. One 
minimizer, u + , is close to u — 1 except for the boundary 
layers at y = and y — 1 (case (i)) or at the edge of 
square (case (ii)). The other minimizer, is close to 
u = — 1 except for the boundary layers at x — and 
x = 1 (case (i)) or at the edge of square(case (ii)). In the 
absence of noise, both minimizers are stable equilibrium 
states. 

In Figure ||, we show the time sequences of the switch- 
ing process for different values of 8 at fixed T. The over- 
all trend is consistent with what was found in the previ- 
ous example, namely there are more and more nucleation 
events as 8, T — * 0. 

Finally let us consider thermally activated switching 
of a magnetic thin film, modeled by the Landau-Lifshitz- 
Gilbert equation, which after suitable nondimcnsional- 
ization, reads 

m t — f(m) := —m x h e jj— am x [m x h e ff). (10) 

Here m is the magnetization distribution (3d vector) 
which satisfies |m| = 1, h e jj is the effective local field, 
given in terms of the free energy by h e g= —8E[m]/ '8m, 
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FIG. 5: Minima and saddle points configurations crossed by 
the optimal switching path from (ml) to (m4); the energy of 
the magnetization during the switching shows that two inter- 
mediate minima, (m2) and (m3), and three saddle configura- 
tions, (si), (s2), and (s3), are crossed during the switching. 
The out-of-plane component of the magnetization vector m is 
very small (less than 10~ 2 )during the switching and we only 
plot its in-plane component: it points to the right in the white 
areas, to the top in the gray areas, and to the left in the black 
areas. 

where 

E[m}= I |Vm| 2 d 3 a;+ / <j>(m)d 3 x + [ \Vu\ 2 d 3 x. 
Jn Jn Jr 3 

Here the three terms represent respectively energies due 
to exchange, anisotropy, and stray field. Q is the region 
occupied by the magnetic thin film. The potential u, 
defined everywhere in space, solves div(— Vw + m) = 0, 
where m is extended as outside the sample. 

The first term at the right hand side of ( [l~0| ) describes 
precession of m around h e jj; the second term, which can 
be written as —am x (m x h e jj) — a(h e jj— (m • h e jj)m), 
is a damping term whose strength is measured by the pa- 
rameter a. We choose physical parameters corresponding 
to permalloy, and consider a thin square domain of size 
200nm x 200nm x lOnm. 

The dynamics in ( jlpj ) has equilibrium states corre- 
sponding to the minima of the magnetic energy, E[m]. 
Thermal noise effects can be represented by modifying 
h e g as h e jj+ ^/e 77, where rj is a vector white-noise, and 
will eventually switch the system from one minimum to 
any other one. In Fig. |^, (ml) and (m4) show two min- 
ima where the magnetization points mostly to the right 
(white color) or to the left (black color), respectively. For 



small noise, the optimal path for a switching from (ml) 
to (m4) during the time interval [0, T] is obtained by 
minimizing the action functional corresponding to ([!(]) : 

S T [m}= [ [ \m t - f(m)\ 2 d 3 xdt. (11) 
Jo Jn 

The global minimizer of ( |i~l"| ) on a long time interval, 
T = 450, is presented in Fig. |5|. Looking at the energy 
E[m] during the switching one sees that two intermediate 
minima, (m2) and (m3), and three saddle configurations, 
(si), (s2), and (s3), are crossed during the switching. The 
first step of the switching from (ml) to (m2) through (si) 
rotates the interior magnetization by 90°. The second 
step from (m2) to (m3) through (s2), which is the most 
expensive step, switches magnetization at the top and 
bottom end domains. Finally, the third step from (m3) 
to (m4) through (s3) is similar to the first one. 

In conclusion, we have shown that the least action prin- 
ciple can be turned into an efficient numerical procedure 
for finding the optimal dynamic path in spatially- ex- 
tended systems driven by small noise, and we have pre- 
sented applications to bistable systems modeled by the 
Ginzburg-Landau equation and the magnetic thin films. 
Our numerical method circumvents the difficulty of hav- 
ing to compute too many time steps in order to observe 
the relevant events. It is general and can be applied to a 
variety of problems in physics, chemistry as well as biol- 
ogy. 
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